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^ : Abstract 

\^ ' We study a model with Major ana neutrino masses generated through doubly charged scalars at 

two-loop level. We give explicit relationships between the neutrino masses and the same sign dilep- 
^^\ ton decays of the doubly charged scalars. In particular, we demonstrate that in the tribimaximal 

D ' limit of the neutrino mixings, the absolute neutrino masses and Majorana phases can be extracted 



through the measurements of the dilepton modes at colliders. 
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I. INTRODUCTION 

It has been revealed that at least two light neutrinos have nonzero masses and the mixing 
matrix is characterized by two large mixing angles from neutrino oscillation experiments. 
These evidences exhibit new physics beyond the standard model (SM). However, the origin 
of neutrino masses is still mysterious even though considerable efforts have been put in 
both theory and experiment for decades. In particular, two crucial properties of neutrinos, 
which can not be disclosed in the neutrino oscillation experiments, are the Majorana nature 
of neutrinos and the absolute values of the neutrino masses. Presently, the only neutrino 
experiment that could provide a direct evidence of Majorana neutrinos is the search for 
the neutrinoless double beta (Oi//3/3) decays, which could be used to set some limits on the 
absolute neutrino masses l|. The absolute neutrino mass can be also measured through 
the electron energy spectrum away from the end-point in the nuclear beta decay (i.e. the 



tritium decay, m/? = ^/^~\Uei^rn^^ < 2 eV) [2|, y] , while the sum of the neutrino masses 
has been constrained from cosmology, given by J2i^^Vi < 0-58 eV (95% CL) [4|. 

Instead of probing neutrino properties at low energy experiments, the possibility to ex- 
plore them at the Large Hadron Collider (LHC) has been widely studied in the literature. 
There are three different types of seesaw mechanisms to generate neutrino masses at tree 
level. The "minimal" type-I seesaw mechanism [5] can not be directly tested at collider ex- 
periments due to the suppressed mixings between the light neutrinos and heavy right-handed 
singlets unless some symmetry is introduced [6|-|9[. However, the unsuppressed gauge interac 



tions [l2|,[l3|of 



3oth scalar and fermionic triplets will help us to test the ideas of type-II [1 



and type-Ill [ll| seesaw mechanisms at the LHC, respectively. Furthermore, an one to one 
correspondence between dilepton decay widths and neutrino masses exists in the type-II 
seesaw model where the triplet scalar T couples to the left-handed lepton doublets El via 
the gauge invariant Yukawa interaction: 

L = taLhabir2T(.bL + h.c. , (1) 

where the indices a, h denote e, /i, r and Ejl represents the left-handed lepton doublet of the 
jth flavor. The neutrino mass matrix is given by 

rriyah = V2habVT (2) 

after the triplet receives the vacuum expectation value (VEV), (T) = vt/\/2. The decay 



width of each same sign dilepton decay mode of the doubly charged Higgs scalar P^^, 
T{P^^ — )■ iaL^b^b)^ is directly related to the corresponding neutrino mass matrix element 
through Eq. (j2]). Thus, the discovery of P^"^ may help us to understand the Majorana 
nature of neutrinos, and by studying the branching ratios of the dilepton channels we may 
obtain some important informations such as the absolute neutrino masses and Majojrana 



phases. The related phenomenologies have been extensively studied in the literature 14 1. 



In this work, we study a model originally proposed in Ref. 15|, in which a discrete 



symmetry is imposed to forbid the Yukawa coupling in Eq. ([T]) at tree level. As a result, the 



neutrino masses are generated at two-loo p le vel with the normal hierarchy spectrum 16| . 



while the OufSP decay arises at tree level 15|, [17(1 . One of the interesting properties of the 



model is that it shares the same feature of the direct link between the decay widths of the 
same sign dilepton modes and neutrino masses as the type-II seesaw model. We will focus 
on the parameters which can not be measured in the neutrino oscillation experiments such 
as the absolute masses of three light neutrinos and the Majorana phases ipi and 1^2- By 



assuming the tribimaximal mixings [18|] and utilizing the measured mass square differences 
for the neutrinos, we derive some explicit relations between the neutrino masses and the 
branching fractions of the doubly charged scalars to the charged lepton pairs. In addition, 
as the neutrino masses are proportional to the products of the charged lepton masses maTrih 
due to the loop integral, the fractions l^\p±±^l±''^±] = ^ and br\p^^Zt^t4 = ^ ^^® 
much larger than those of unity predicted in the type-II seesaw model in the limit of the 
tribimaximal mixings. Clearly, we are able to differentiate the two types of the models by 
counting the events arising from the dilepton decays of the doubly charged Higgs scalars at 
the LHC. 

Our paper is organized as follows. In Sec. II, we briefly introduce the model. In Sec. Ill, 
we relate the neutrino masses with the dilepton modes of the doubly charged scalars. We 
conclude our results in Sec. IV. 



II. THE MODEL 

The simplest version of the model has been given in Ref. fl5l with some of its phe- 

n 

nomenologies presented in Ref. [16|. The idea of the model is to suppress the Yukawa 
interaction in Eq. ([1]) at tree level and induce it radiatively. Here, we give an explicit exam- 



plqj to forbid the tree contributions to the neutrino masses. The model consists two Higgs 
doublets 01 and 02, one triplet T and one doubly charged singlet "^ with the hypercharges of 
—1/2, —1 and 2, respectively. Besides the gauge symmetry, there is a Z2 discrete symmetry 
with the transformations of (pi — ?■ +0i, 02 ^ —02, T — )■ — T, ^ — )■ +\1', and f -^ f, where 
/ represents the SM fermion. The most general potential in this model can be written as 

V = -/X?0l0l + Ai(0l0i)2 - fil4(P, + A2(0^02)' 

-l4Tr{T^T) + XT[Tr{T^T)]'^ + \'j,Tr{T^TT^T) 

+K^,Tr{4>\<p{T^T) + «:;^01ttVi + «:„,0l0i^t^ 

+«:^,Tr(0|02rtr) + «:;^0^TTt02 + «:,,,0^02^t* 

+A30l0i0|02 + A40l020^0i + pTr{T^Tm^^) 

+ (M0[tV2 + A50I 020102 + A0lT0;^ + H.C.), (3) 

where 0j = «T20*. Under the symmetries, the doubly charged singlet \1' couples to the 
right-handed charged leptons Ir via the Yukawa interaction, 

Ly = Yabi^adbR^ + h.C, (4) 

where Yah is a 3 x 3 symmetric matrix with the indices a, h stand for e, /i, r and Ijr represent 
the right-handed lepton singlets. We note that the interaction between the scalar triplet and 
the left-handed lepton doublets is not allowed by imposing the Z2 symmetry such that there 
is no neutrino mass term at tree level unlike the type-II seesaw mechanism. When the scalar 
fields 01^2 and T develop VEVs, both the gauge and Z2 discrete symmetries spontaneously 
break down. The neutrino masses will be generated through two-loop diagrams as shown in 
Ref. [l5j, given by 

[mujab = /= maTTlh [I{MpJ - /(MpJJ , (5) 

where w^ < 4 GeV is the VEV of the scalar triplet T, bounded by the p-parameter (= 
My^/Mz cos^ 9w) [2, ll9[, rna^b stand for the charged lepton masses, Mp^,^ are the masses of 
the doubly charged scalar eigenstates Pi 2 with the mixing angle 6 defined by 

P-t^ I \ -sin^ cos^ / I ^±± 



^ Two possible scenarios have been mentioned in the footnote of Ref. 



i5i. 



and the integral I{Mp) is expressed as 

liM- )= [^ [ -^^ ^ ^ 1 1 (7) 



which can be approximated to be J ~ u^dj^^i log^ I j^ ) for Mp. > M^ [20|- From Eq. 
(|5]), we see that the one to one correspondence between neutrino mass elements m,^^^ and the 
Yukawa couplings Yab provides us the opportunity to determine the neutrino masses through 
the measurements of the dilepton modes. We note that even for small neutrino masses, the 
OuPP decays in this model can be large since they are dominated by the exchanges of the 



doubly charged scalar s at tree level |l5 |. 



III. DECAY BRANCHING RATIOS RELATED TO NEUTRINO MASSES 

A. Dilepton decays 

There are several channels that the doubly charged scalars can decay into, such as P^^ — >■ 
^aR^bw W^W^, W^P^, P^P^ and W^W^P^, where P±± are referred to the lighter mass 
eigenstate among P^^^ {i = 1,2). Due to the kinematical consideration, we expect that 
the first two kinds of the modes contribute the most part to the width of P^^ if Mp±± ~ 
Mp± ~ Mpo. However, as long as the mass splitting between P^^ and P^(Mpo) is large 
enough, the last three channels may dominate at the high mass region. In our discussion, 
we assume that the last three types of decay channels with scalar(s) in the final states are 
not allowed. Hence, the two-loop suppression factor in neutrino masses of Eq. (J5]) makes 
the fraction of ^jflty ~ ^^^ ( '"'m.C'^' )^ i^X ~ ^'(lO^O), so that the branching ratio 
with the dilepton final states is almost 100% for vt < 0{1) GeV by taking m.y^^ = 0.1 eV 
and Mp^ = 200 GeV. While the fraction in the seesaw type-II model is {^^f{^)^ > 1 
only for vt < C(10~^) GeV with the same input. We take the values of m,^^^ ~ 0.01 eV, 
"^i^eM ~ ^^feT ~ 0.1 eV, and rriu^^ ~ "^i/^r ~ "^i^rr ~ 1 eV to illustrate the possible branching 
ratios of the dilepton modes in our model. These values are based on the texture of the 



neutrino mass matrix with the normal hierarchical spectrum as predicted by this model 16 1. 
The estimations of the decays for Mp = 200 GeV and f ^ = 1 GeV are shown in Table [B 
One should keep in mind that there might exist some cancellations among the combinations 
of three neutrino masses and Majorana phases in each element of the neutrino mass matrix. 



For example, the component rriy^^ may go to zero with certain values of Majoran phases and 
the lightest neutrino mass for the normal hierarchical spectrum. In this case, Oz//3/3 decays 
may be out of the reach of the current experimental sensitivity. Therefore, the branching 
ratios of the dilepton modes shown in Table [I] will change drastically with different values 
of Majorana phases and neutrino masses, which will be discussed in Sec. Ill C. 

The production of P"^^ is dominated by the Drell-Yan process, for which the next to 
leading-order contribution from QCD enhances the cross section by a factor of 1.25 at the 
LHC 21j. Several simulations have been performed for P^^ — )■ i^^i^ at the LHC 22j. For 
channels of the doubly charged Higgs decays into e or fi, the SM background is shown to 
be negligible in the signal region of the high invariant mass close to Mp±±. In contrast, the 
spectrum of the missing transverse momentum for the r decay product will be softer due 
to the subsequent decays of r — )■ eui? and r — )■ fiuu with a branching ratio around 17% in 
each channel. Additional backgrounds with jets, such as W^W^jj, may fake the hadronic 
T decays. As a result, the r tag efficiency is about 50% and the fake rate is around 1% |23|. 
The detector-specific error analysis is beyond the scope of this paper. We have used uniform 
uncertainties for all branching ratios for the rough estimation of the effect. For example, the 
approximate event numbers for the pair production of P^^ are 15000, 3000 and 900 with 
Mp = 200, 300 and 400 GeV, respectively, at the LHC (v^ = 14 TeV) for the integrated 
luminosity of L = 300fb-\ BR{P^^ -^ i^i^) = 100%, and the efficiency e^jf = 0.5. We 
write the dilepton decay widths as 

IV |2 



r(p 



~^ ^aR^bW 



^7r(l + 6ab) 



siMp. 



(8) 



TABLE I. The branching ratios of P — t- ^^^f, i^ab ~ ^ '/^ ^^"^ ''" ) by assuming the neutrino 
mass elements of rriv^^ = 0.01, ruy^^ = m^^^ =0.1, and ruu^^ = rriy^^ = niy^^ = 1 eV for vt = I 
GeV and Mp = 200 GeV []j. 



BRe±e± 


BRe±^± 


BR^±^± 


^^^±M± 


BRf,±r± 


BR^±^± 


0.995 


4.6 X 10-3 


1.6 X 10"^ 


5.4 X 10"^ 


3.8 X 10"^ 


6.7 X IQ-ii 



B. Neutrino masses and mixings 



Since the decay widths in Eq. (|8]) are proportional to Yat, which are related to the neutrino 
masses in Eq. i^ it provides with us an opportunity to study the spectrum of neutrinos at 
the LHC We can express the branching ratios of the same sign charged lepton pair modes 
in terms of the components of the neutrino mass matrix 



BR 



■ab 



^i^B^. 



bRJ 



siug Mpi 



X 



\7n 



Vab\ 



total 



Ttotai X A-ng^vl sin^ 29 [l{Ml^ - I{Ml^))\ (1 + SabWafnl 



\7n 



oc ^—^ (9) 

It will be clear later that the overall factor including the total decay width is irrelevant and 
the dependence of charged lepton masses appearing in the dilepton branching ratios is due 
to the loop integral in our model. 

Similar to the CKM mixing matrix in the quark sector, the neutrino mass matrix can be 
diagonalized by the unitary matrix Upmns, defined by |24| . 

'mi ^ 

mae^^i 
mse'"^^ 

where ?/'i_2 are referred to as the Majorana phases and the PMNS matrix can be parametrized 

as 



M„ = U, 



PMNS 



\ 



^ PMNS^ 



(10) 
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PMNS 



C12C13 



S12C13 



\ 



-S12C23 - Ci2S23Sl3e*'' C12C23 - Si2S23Si3e*'' 



v 



A& 



AS 



(11) 



/ 



si3e 
S23C13 

S12S23 - 0120235136'" -C12S23 - 5120235136"" C23C13 

with the Dirac phase of 5. Currently, the neutrino oscillation experiments can only give the 
mass-square differences |2[: 

Am^^ = (7.59 ± 0.20) x 10"^eV^ jAm^al = (2-43 ± 0.13) x 10"^eV^ 



(12) 



and the mixing angles [2|: 



sin^ (2^12) = 0.87 ± 0.03, sin^ (2^23) ^ 1, sin^ (2^13) < 0.19. 



(13) 



A very good approximation of leptonic mixing matrix is proposed with the so-called tribi 
maximal mixing form 18 1 

i sin' ^23 =2 



sin^e 



12 



sin^^ 



13 



0. 



(14) 



In this limit, we can express the elements of the neutrino mass matrix as 

4 4 1 

I'^i-eeP = q"^! + q"ii"^2COS?/^i + -ml , 

12 1 

l^^i^ej^ = g"^i - g"^l"^2Cos?/'l + -ml, 

I |2 I |2 

I i^er I I '^e/j, I " 

I'^'^MM 1^ = ^"^1 + n"^2 + 7"^3 + n"^l"^2 COS t/^i + -mim^ COS 7/'2 + -m2m3 COS (V'l - ^2) 

00 y 4 y o 

l"^i^urP = 7T7:"^i + 7t"^2 + 7"^3 + -m,im2 cos-i/'i - -mim3COS^2 - 7rm2m3 cos (■?/'i -V'2) 

do y 4 y 3 



muj'. (15) 



|2 _ |_ |2 



C. Relations 

Since the normal hierarchical mass spectrum of the lisrht neutrinos is predicted in our 

n 

model pj(, we can parametrize the eigenvalues of the masses in terms of the mass differences 
and the lightest one m,i, given by 



mi, m,2 = \/ m,1 + Amli, 7713 = Wmf + Am^i + Am|2- (16) 



By utilizing Eqs. (19]) and ( jTSl) . we can define the quantity 

ml{2mlBR^^ + mlBR^, + mjBR,^) ^ ml + fAmJi + |Am^, 
^ 2ml{mlBRee + mlBR^^,) m? + iAm^^ ^ -^ 

such that we can determine the lightest neutrino mass mi via the quantity Ci by measuring 
the branching ratios of the ee, e/z, /x/i and /xr channels with the relation 

(I - \Ci)/\ml, + lAm^, 
mi - ^r^l • (1«) 

Similarly, we can express the masses m2 and m^ as 

(C2 + 2)/\ml^ + SCaAm^ 
3(1 - C2) 

and 



m^ = 1-2 ^ ^;^»^2i ^ o^2^»^32 (ig) 



(9 - 8C3)Am|2 + (6 - C3)Am2 
9 - 17^3 



^2 ^ w-"^3;^"^32^l^-^3;^"^2i^ (20) 



with the quantities C2 and C3, defined by 






8 



and 

' ml{mlBRee + mlBR,^) ^ ' 

respectively. As a result, we are able to indirectly determine the neutrino masses by measur- 
ing the quantities Ci, C2, and C3 in the tribimaximal limit of the neutrino mixings. One can 
see that the different terms in numerators of Cj(i =1-3) are all of the same order of magni- 
tude since the large differences in the branching ratios (see Table [T]) are compensated by the 
large differences of the masses of the charged leptons. The degree of accuracy in measuring 
Ci depends on the assumption of a sufficient number of like-sign leptons to be observed, so 
the sensitivity depends on the number of P^"^ and the branching ratios of dilepton modes. 
As we mentioned in Sec. Ill A, the branching ratios of the dilepton channels are sensitive 
to the values of ipi 2 and mi as given in Eq. ( TTSl) . Thus, BR,±,± will be very different from 
those shown in Table |T] with different values of ipi^2 and nii. However, we would use the ratio 



BRee lmf^\m„^J'^ 1 m^ |mf + fmi^m^ + Am^^ cost/'i + |Am: 



BRe^, 2 ml \m^,X 2 ml |m? - |mi v^mf + Am|i cos V^i + ^A 



(23) 



to pin down the allowed parameter region as a function of the lightest neutrino mass m\. 
This is demonstrated in Fig. [T] (right). Similar results of other ratios, such as ■§^^, ( Jr^) 
||2i^), and (||2^, If^) are displayed in Figs. [D (left), g] (left, right), and E] (left, right), 
respectively. It is interesting to note that the branching ratio of P^^ — )■ e^e^ reduces 
significantly in the small region around mi ~ 0.005 eV due to the cancellations. In this 
region, the branching ratios of the rest dilepton modes will be enhanced (see Eq. fl23|) . 
Figs. [H and Fig. [2]). Our relations of Eqs. (TT7|) -( 122|) will be practically implemented for the 
lowest mass of the doubly charged Higgs scalars. If there is no cancellation in viiy^^ , we can 
still try to narrow down the parameter space of mix and 't\}\^2 by measuring the fractions of 
the dilepton modes. For instance, in the limit of mi\ — )■ the ratio of BR^e/ BR^^ becomes 

BRe^ = ^. (24) 

BR,^ 2ml 

On the other hand, if mii is measured, the Majorana phase ipi can be extracted from the 
relation 

{7C,-?,)m\ + 2C,/\ml, 
cosipi = — — — — , (25) 
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FIG. 1. Ratios of BR^e/BR^^ (left) and BRee/BR^fj, (right) versus the hghtest neutrino mass rrii 
with scanning over the possible values of Majorana phases, where the shadow areas are the allowed 
regions. 
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FIG. 2. Legend is the same as Fig. [T] but with the ratios of BR^e/BRer (left) and BR^^/BR^^ 
(right). 



where C4 is expressed as 



C4 



2mlBRee - mlBRe^ 
2mlBRee + mlBRe 



(26) 



In Fig. m we plot the allowed region of ipi versus 



i3i\e.e — i^ tie 



D„ , „„ and mi. Note that the phase ihi 
becomes indefinite in Eq. f l25|) if we set mi — ?■ 0. This is because when the lightest neutrino 
mass is zero, there is only one Majorana phase left, related to the relative phase of ipi — 1^2 
as shown in Eq. ( fTSl) . In this limit, we obtain 



2/BRu 



COS (V^i - ij2) 



3 ^ BRe 



mi ^i 



i|Am2 + |Am2 



A?7i2i vAm|7+Am 



(27) 



32 



The allowed region of ipi — ip2 for mi — )■ is displayed in Fig. [5l where we havej;aken the 
uncertainties of the mass differences from the the solar and atmospherical data 
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FIG. 3. Legend is the same as Fig. [T]but with the ratios of BRe^/BR^r (left) and BR^^^/BR^r 
(right). 



ri(eV) 




FIG. 4. The Majorana phase ipi in terms of C4 = 9^^rb''^_i_^!^rb'''" and mi. 
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IV. CONCLUSION 

We have studied the close relationships between the neutrino masses and the same sign 
dilepton decays of the doubly charged Higgs scalars in the model in which the neutrinos 
are Majorana particles with their masses are generated radiatively at two-loop level. Since 
the dilepton modes in our model could be reachable at the LHC, it is natural to use their 
branching ratios to infer the neutrino masses. We have explicitly shown that in the limit of 
the tribimaximal mixings, the absolute scale of neutrino masses can be expressed in terms of 
Ci (i = 1-3) based on the certain combinations of dilepton branching ratios. It is possible to 
determine the neutrino masses by just counting the events arising from the dilepton decays 
of the doubly charged Higgs in its lowest mass region. The allowed parameter space of the 
fractions among each dilepton branching ratio as the function of the lightest neutrino mass 
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FIG. 5. The relative Majorana phase ■01 — 02 versus „„"'' for m-i — )• 



mi is presented. These relations combined with the data from other neutrino experiments 
may help to set a limit of mi and Majorana phase 01^2 in the future. 

Acknowledgments 

This work is supported in part by the National Science Council of ROC under Grant 
#s: NSC97-2112-M-006-004-MY3 (CSC), NSC-95-2112-M-007-059-MY3 (CQG) and NSC- 
98-2112-M-007-008-MY3 (CQG) and by the Boost Program of NTHU (CQG). 



[1] S. M. Bilenky, S. Pascoli and S. T. Petcov, Phys. Rev. D 64, 053010 (2001) 
|arXiv:hep-ph701022 65 1 ; S. Pascoli, S. T. Petcov and T. Schwetz, Nucl. Phys. B 734, 24 
(2006) | arXiv:hep-ph/0505226] . 

[2] C. Amsler et al. (Particle Data Group), Physics Letters B667, 1 (2008) and 2009 partial update 
for the 2010 edition. 

[3] E. W. Otten and C. Weinheimer, Rept. Prog. Phys. 71 (2008) 086201 |arXiv:0909.2T04] [heo- 
ex]]. 

[4] E. Komatsu et al., |ar2£iYilOOL4538l [astro-ph.CO]. 

[5] P. Minkowski, Phys. Lett. B 67, 421 (1977); M. Gell-Mann, P. Ramond, and R Slansky, Pro- 
ceedings of the Supergravity Stony Brook Workshop, New York, 1979; T. Yanagida, Proceed- 
ings of the Workshop on the Baryon Number of the Universe and Unified Theories, Tsukuba, 
Japan, 1979; R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

12 



[6] A. Datta, M. Guchait and A. Pilaftsis, Phys. Rev. D 50, 3195 (1994) [arXiv:hep-ph/9311257] . 
[7] F. M. L. Almeida, Y. D. A. Coutinho, J. A. Martins Simoes and M. A. B. do Vale, Phys. Rev. 



D 62, 075004 (2000) larXiv:hep-ph/00020241. 



[8] T. Han and B. Zhang, Phys. Rev. Lett. 97, 171804 (2006) [arXiv:hep-ph/0604064 1 ; P. Fileviez 

Perez, T. Han and T. Li, Phys. Rev. D 80, 073015 (2009) [arXiv:0907.4T86] [hep-oh]]. 
[9] F. del Aguila, J. A. Aguilar-Saavedra and R. Pittau, JHEP 0710, 047 (2007) 

|arXiv:hep-ph/0703261| ; F. del Aguila and J. A. Aguilar-Saavedra, Nucl. Phys. B 813, 22 

(2009) [ arXiv: 0808. 24681 [hep-ph]]. 
[10] M. Magg and C. Wetterich, Phys. Lett. B 94, 61 (1980); J. Schechter and J. W. F. Valle, 

Phys. Rev. D 22, 2227 (1980); T. P. Cheng and L. F. Li, Phys. Rev. D 22, 2860 (1980); 

G. B. Gelmini and M. Roncadelh, Phys. Lett. B 99, 411 (1981); G. Lazarides, Q. Shah and 

C. Wetterich, Nucl. Phys. B 181, 287 (1981); R. N. Mohapatra and G. Senjanovic, Phys. Rev. 

D 23, 165 (1981); J. Schechter and J. W. F. Valle, Phys. Rev. D 25, 774 (1982). 
[11] R. Foot, H. Lew, X. G. He and G. C. Joshi, Z. Phys. C 44, 441 (1989); E. Ma, Phys. Rev. 



Lett. 81, 1171 (1998) (arXiv:hep-ph/9805 219|. 

[12] K. Huitu, J. Maalampi, A. Pietila and M. Raidal, Nucl. Phys. B 487, 27 (1997) 
|arXiv:hep-ph/9606311] ; A. G. Akeroyd and M. Aoki, Phys. Rev. D 72, 035011 (2005) 
|arXiv:hep-ph/0506176] ; P. Fileviez Perez, T. Han, G. Y. Huang, T. Li and K. Wang, Phys. 
Rev. D 78, 071301 (2008) ta rXiv:0803.3450] [hep-ph]]: P. Fileviez Perez, T. Han, G. y. Huang, 
T. Li and K. Wang, Phys. Rev. D 78, 015018 (2008) [arXiv:0805.3536 [hep-ph]]. 

[13] R. Franceschini, T. Hambye and A. Strumia, Phys. Rev. D 78, 033002 (2008) |arXiv:0805. 16131 
[hep-ph]]. 

[14] J. Garayoa and T. Schwetz, JHEP 0803, 009 (2008) |arXiv:0712.1453. [hep-ph]]; 
A. G. Akeroyd, M. Aoki and H. Sugiyama, Phys. Rev. D 77, 075010 (2008) [ arXiv:0712.40T9] 
[hep-ph]]; M. Kadastik, M. Raidal and L. Rebane, Phys. Rev. D 77, 115023 (2008) 
|arXiv:0712.3"9T2] [hep-ph]]; C. S. Chen, C. Q. Geng and D. V. Zhuridov, Phys. Lett. B 666, 
340 (2008) [a rXiv:0801.20TT] [hep-ph]]: W. Chao, Z. G. Si, Z. z. Xing and S. Zhou, Phys. Lett. 



B 666, 451 (2008) [arXiv:0804.1265 [hep-ph]]; P. Ren and Z. z. Xing, Phys. Lett. B 666, 48 

(2008) [ arXiv:0805.4292] [hep-ph]]: H. Nishiura and T. Fukuyama, Phys. Rev. D 80, 017302 

(2009) |arXiv:0905.3963] [hep-ph]]. 

[15] C. S. Chen, C. Q. Geng and J. N. Ng, Phys. Rev. D 75, 053004 (2007) |arXiv:hep-ph/0610118| 

13 



[16] C. S. Chen, C. Q. Geng, J. N. Ng and J. M. S. Wu, JHEP 0708, 022 (2007) [arXiv: 0706 .1964] 

[hep-ph]]. 
[17] J. Schechter and J. W. F. Valle, Phys. Rev. D 25, 2951 (1982). 
[18] P. F. Harrison, D. H. Perkins and W. G. Scott, Phys. Lett. B 530, 167 (2002) 

|arXiv:hep-ph/0202074] ; X. G. He and A. Zee, Phys. Lett. B 560, 87 (2003) 

|arXiv:hep-ph/0301092] . 
[19] J. F. Gunion, R. Vega and J. Wudka, Phys. Rev. D 43, 2322 (1991). 
[20] J. van der Bij and M. J. G. Veltman, Nucl. Phys. B231, 205 (1984); K. S. Babu and 

C. Macesanu, Phys. Rev. D 67, 073010 (2003); K. L. McDonald and B. H. J. McKellar, 



arXiv:hep-ph/0309270 



[21] M. Muhlleitner and M. Spira, Phys. Rev. D 68, 117701 (2003) [arXiv:hep-ph/0305288] . 



[22] G. Azuelos, K. Benslama and J. Ferland, J. Phys. G 32, 73 (2006) l arXiv:hep- ph/0503096 1 ; 
T. Han, B. Mukhopadhyaya, Z. Si and K. Wang, Phys. Rev. D 76, 075013 (2007) 
|arXiv:0706i0441 [hep-ph]]; A. Hektor, M. Kadastik, M. Muntel, M. Raidal and L. Re- 
bane, Nucl. Phys. B 787, 198 (2007) |arXiv:0705. 14951 [hep-ph]]; B. C. AUanach et al, 



arXiv:hep-ph/0602 198 [ Part 28 (A. Hektor et al). 



[23] CMS TBR: CMS physics: Technical Design Report V.2 : Physics Performance, CERN-LHCC- 
2006-021; ATLAS TDR: ATLAS detector and physics performance. Technical design report. 
Vol.2, CERN-LHCC-99-15. 

[24] Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 28, 870 (1962); B. Pontecorvo, Sov. 
Phys. JETP 26, 984 (1968) [Zh. Eksp. Teor. Fiz. 53, 1717 (1967)]. 



14 



